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Abstract 

Several inequalities are proved for the mixing time of discrete- 
time quantum walks on finite graphs. The mixing time is de- 
fined differently than in Aharonov, Ambainis, Kempe and Vazirani 
(2001) and it is found that for particular examples of walks on a 
cycle, a hypercube and a complete graph, quantum walks provide 
no speed-up in mixing over the classical counterparts. In addi- 
tion, non-unitary quantum walks (i.e., walks with decoherence) 
are considered and a criterion for their convergence to the unique 
stationary distribution is derived. 

1. Introduction 

The origin of the concept of quantum walk lies in quantum computation theory, 
where a quantum version of the classical random walk was invented in an attempt to 
improve over classical computational algorithms. The early papers that formulated 
the main ideas of quantum walk are and ifTTTl . Among numerous later papers, we 
would like to point out Q where the continuous-time quantum walk was defined 
and HI which defined and studied the discrete-time quantum walk on finite graphs. 
An introductory review of quantum walks written from the prospective of quantum 
computation can be found in Q. For recent developments the reader can also 
consult GO. 

From the beginning, it became clear that quantum walks on both finite and infi- 
nite graphs have many differences from the classical walk. For example, the prob- 
ability to find a particle at a particular vertex of a finite graph does not converge to 
a limit but in general oscillates forever. However, the average of this probability 
over time does converge to a limit, which can be interpreted as follows. We start 
quantum walk in a certain state and measure the particle at a random time t, which 
is distributed uniformly over interval [0, T] . This measurement finds the particle 
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at a particular vertex v with a probability p(v,T) , which converges to a limit as 
T — > co. How large should T be if we want to make sure that p (v, T) is close to 
its limit? 

Let us introduce some definitions to make this question more precise. 

The quantum walk on a finite graph is a 4-tuple (G, S,ip,U) , where G = ( V, E) 
is a finite graph, 5 is a finite set, ip is a function in (V x S), and [/ is a unitary 
operator on L 2 (V x S) . It is assumed that [|^[| = 1. Elements of S are called 
chiralities and the function ip t = is the wave function at time t G Z. If a 
measurement is performed over the system at time t, then the walking particle is 
found at vertex v in state s with probability \ip t (v, s)\ 2 . 

We assume that the quantum walk is local. That is, let x and x' denote pairs 
(v, s) and (V, s') , respectively. A quantum walk is local if U x > x = (5 x i, U5 X ) ^ 
implies that v ~ v', that is, vertices v and v' are connected to each other. A local 
quantum walk is called the general quantum walk in U. 

A special case of the general quantum walk is the coined quantum walk (ifTI). 
Here is how it is defined. Let G be a d-regular graph and let S = {1, . . . , d} . 
Assume that the neighbors of each vertex v are labelled as vi where i = l,...,d. 
In addition, assume that if v ^ w and V{ = Wj then i / j. (Such a labelling L 
always exists on Cayley graphs of finitely-generated groups, where we can identify 
elements of S with generators and inverses of generators of the group and write 
v 9 = vg and v g -i = vg^ 1 . In this case the choice of labelling is equivalent to the 
choice of ordering of generators and their inverses.) 

Define U as follows. Let x and x' denote pairs (v, s) and (v', s') , respectively. 
If v' <yo v, then U x i x = 0. Otherwise, v' = vi and U x > x = 5i s C s i s , where C is a 
unitary matrix acting on I? (S) , which is called the coin of the quantum walk. It 
is easy to check that matrix U is unitary. Intuitively, let the particle be at vertex 
v in state |s) . Then, at the next moment the particle will be at vertex v s in the 
superposition state C \s) . This is the coined quantum walk on G corresponding to 
labelling L and coin C. 

A typical example of the coined quantum walk is the Hadamard quantum walk 
on the cycle Z n . In this case, the coin is the Hadamard transformation: 




Another popular choice of the coin is Graver's transformation: 
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That is, the state s remains unchanged with amplitude (l — |) and moves to s' ^ s 
with amplitude —2/d. We will call walks with this coin the Grover quantum walks. 

A generalization of this concept is the non-unitary quantum walk [1|. A non- 
unitary quantum walk is specified by 4-tuple (G, S, p, T) , where G and S are as 
before a finite graph and a finite set, p is a density matrix (i.e., a positive unit-trace 
operator on (V x 5)), and Tis a completely-positive trace preserving operator 
acting on density matrices. In the literature, T is called a superoperator lfl4ll . or 
a quantum channel lfl4l . or a trace-preserving quantum operation HI 31 , We will 
use these terms as synonyms. Let x denote a pair (v, s). The probability to find a 
particle at vertex v in state s at time t is given by (x\Tp\x) . A non-unitary quantum 
walk is local if (x'\T (p) \x') > for x' = (V, s') implies that there is x = (f , s) 
with v ~ v' such that > 0. (The concept of locality is more complicated in 

the non-unitary case and this definition is different from the definition in iTTTl . ~> 

An example of a non-unitary quantum walk is given by a weighted sum of uni- 
tary quantum walks. In this example, T (p) = Yli=iPiUipU* where U{ are uni- 
tary operators, pi > 0, and Y2i=iPi = 1- Intuitively, an operator Ui is used at each 
step of the walk with probability pi . If all U{ are local, then T is also local. Another 
example is T (p) = pY^\=i PipPi + {I — p)U pU*, where Pi are projections and 
Yli=i Pi = P This is a walk in which with probability p the particle is measured 
and with probability 1 — p it is evolved according to the unitary operator U. 

First, let us consider the case of unitary quantum walks. The probability dis- 
tribution \ip t (a?) | 2 in general does not converge to any particular limit. Indeed, all 
eigenvalues of the matrix U have unit absolute value. As a consequence, every 
eigenvector of U corresponds to a stationary probability distribution. If the initial 
wave function if) is a non-trivial supeiposition of the eigenvectors with different 
eigenvalues, then ip t continues to oscillate indefinitely. In the classical case this 
phenomenon occurs only when the random walk corresponds to a periodic Markov 
chain, and this case is not typical. 

The time averages of the probabilities \tp t (x)\ 2 do converge, and the limit 

p{x) = lim - V \ifn (x)\ 2 

T— too 1 *■ — ' 

i=0 

exists although may depend on the initial function ip. We will call this limit the 
time-averaged probability distribution of the particle. In order to quantify the con- 
vergence of the initial distribution to this limit, let us define the distance of the 
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initial distribution from its time-average by the formula 

1 T-l 

d(T,4>)= £ 



- ^2 ( x )\ 2 -p( x ) ■ 

i6VxS t=o 

This is the total variation distance between the averaged probability distribution 
at time T and its limit. By analogy with the classical case, the mixing time of a 
general quantum walk is defined as follows: 

t mix ( e ) = sup inf {T : d (T', < e for all T' > T} . 

That is, this is the minimal time which is needed to reduce the distance between 
the worst initial distribution and its time-averaged limit to a quantity less than e. 

Another definition of the mixing time restricts the choice of initial wave func- 
tions. Namely, 

t~ (e) = sup inf {T : d (T') < e for all T' >T}, 

where B is the set of basis states, that is, ip G B if \ip\ 2 is a delta-function con- 
centrated at (v, s) This is the definition used in H). Clearly, t m i x (e) < t m i x (e). 
In the case of classical random walks, these two mixing times are always equal to 
each other. In the the case of quantum walks, they can be different. 

We intend to estimate the mixing time in terms of the distance between eigen- 
values of U. Let Xk = e l ^ k , k = 1, . . . ,m, be the distinct eigenvalues of U. We 
define the distance between Afe and A; as the smallest distance along the unit circle: 

d(X k ,\i) = min{|/3 fc - fa + 2im\ , n £ Z} . 

The relaxation time of operator U is defined as 

t re i = maxd(A fc , A/) -1 . 

k^l 

Finally, let us define the overlap of two functions ip and ip by the formula 

Q( ( p,^)= 0*0 ^ 0*01 • 

xeVxS 

Note that if tp and ip are two wave functions, then < Q < 1 by the Cauchy- 
Schwartz inequality. 

Theorem 1.1. Let U be the unitary transformation on L 2 (V x S) associated with 
a discrete-time quantum walk (not necessarily coined). Let U have m distinct 
eigenvalues and the relaxation time t re \. Then, 

t m ix (e) < 27r log (2m) 
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(The proofs of all theorems are in Appendix.) 

It is interesting to compare this bound with the corresponding result for the clas- 
sical random walk, where t m i x is) < c log ( — - — ) t re i , where 7r m i n is the smallest 
probability in the limit distribution (see for example Theorems 12.3 and 12.4 on 
p. 155 in OH). In many cases the limit distribution is uniform and this bound can 
be written as t m i x (e) < clog (n/e) t re i, where n is the number of vertices in the 
graph. Note, however, that t Te \ have a different meaning in the classical case where 
it denotes the inverse of the difference between 1 (the largest eigenvalue) and the 
second largest eigenvalue (i.e., the inverse of the "spectral gap"). 

Another significant difference in the formulas for the mixing time is that e enters 
as log £ and in the classical and quantum cases, respectively. This is due to the 
fact that the convergence is exponentially fast in the classical case and polynomial 
(even linear) in the quantum case. 

Finally, it is worthwhile to note that in many cases the classical bound t m ix < 
clog (n) t re i is not optimal, and a large literature is devoted to improvement of this 
result to t m ix < drei with a sharp constant c. 

For the lower bound we prove the following result. 

Theorem 1.2. Let U be the unitary transformation on I? (V x S) associated with 
a discrete-time quantum walk (not necessarily coined). Suppose that U has only 
real eigenvectors. Let A and A' be two distinct eigenvalues with the corresponding 
eigenvectors ip and ip'. Assume that d (A, A') < 2 and e < Q (ip, tjj') /80. Then, 

t • (e) > 9M±fl 

In particular, if A and A' are two eigenvalues with the smallest distance between 
them along the circle, then d(X,X') 1 = t re i and we obtain the estimate 

tmix 0) > \Q — 
o £ 

valid for all sufficiently small e. 

The main message of Theorems l 1 . 1 l and fOl is that the relation of the mixing and 
relaxation times in the quantum case is similar to the analogous relation in the clas- 
sical case. However, the relaxation time is defined differently in the quantum case. 
It is not the inverse of the difference between the largest and the second largest 
eigenvalue, but the inverse of the minimal distance between all distinct eigenval- 
ues. 

Previously, the speed of convergence of (unitary) discrete-time quantum walks 
was investigated in HI. The upper bound for the quantum walks that we obtain 
in Theorem 11-11 is similar to the bound in Theorem 6.1 of [1]. The mixing time 
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is O (t re i log (m)) where t re i is the inverse of the minimal distance between the 
distinct eigenvalues of the matrix U. The main difference of our result from the 
result in [H is that we have log (m) instead of log(n), where m and n are the 
numbers of distinct and all eigenvalues, respectively. 

This difference is significant for the case of the discrete walk on the hypercube, 
where the number of eigenvalues is 2 d and the number of distinct eigenvalues is 
d+ 1. In particular, we show that the mixing time on the hypercube is 0(n log n/e) 
and not exponential as was suggested in [ 12] based on previous estimates in (TJ. 

The lower bound that we obtain is in terms of the relaxation time t re i. It essen- 
tially says that the mixing time is Q (t re i) . This bound is different from the bound 
obtained in [ 1 ], which is formulated in terms of a geometrical property of the under- 
lying graph. In addition, the mixing time is defined differently in ifTTl. As a result, 
the mixing time of the Hadamard walk on the cycle is of order O (n log n/e 3 ) in 
Theorem 4.2 of [1], and of order O (n 2 log n/e) in our Example 1. In the classical 
case, the mixing time is of the order O (n 2 log (e -1 )) • 

Now let us consider non-unitary quantum walks. The study of these walks helps 
us to understand how the decoherence affects performance of quantum algorithms. 
It was noted (see (71) that decoherence in quantum walks can be useful for quan- 
tum algorithms. In particular, it appears that a small amount of decoherence can 
speed up the mixing of the walk. Numeric evidence in (7J was later corroborated 
by analytical estimates in lfT5l . More information about decoherence in quantum 
walks and additional references can be found in the review article 0. 

Let M. denote the linear space of Hermitian linear operators acting on (V x S) . 

1 /2 

The space M. is a Hilbert space with respect to the norm ||p|| 2 = [Tr (p 2 )] 
(which we call L 2 -norm). Other useful norms on M are \\p\\ l = Tr (|p|) where 

\p\ = \/~P*' an d IIpIIl 2 (^) = (w 2 )] 1 ^ 2 > where p is a density matrix. We call 
these norms the trace and L 2 (p) norms, respectively. Superoperators are operators 
on Ai which possess some additional properties. Some well-known properties of 
superoperators are summarized in the proposition below. 

Proposition 1.3. Superoperator T is a contraction in the trace norm (i.e., ||7"(p)||i < 
Hpl^J. There exists a density matrix p s t such that Tp s t = Pst- 

This proposition is an immediate consequence of Theorem 9.2 and Exercise 9.9 
in ED. 

Note that in many cases T is not self-adjoint in I? norm. Moreover, recall that 
in the classical case the stochastic matrix of a random walk is always self-adjoint 
with respect to the norm L 2 (p) , where p is the stationary probability distribution. 
(This result can be traced to the fact that every random walk is a reversible Markov 
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chain.) In contrast, the superoperator of a non-unitary quantum walk T is not 
necessarily self-adjoint with respect to the norm L? (p s t)- m fact, it appears that T 
is not even a normal operator (i.e., T*T 7^ TT*) in many situations of interest. 

Proposition 11.31 establishes the existence of the stationary density matrix. How- 
ever, it does not say anything about the uniqueness or convergence properties, and 
we cannot expect that these properties hold in general. For example, a unitary 
quantum walk typically has many stationary density matrices and the convergence 
fails unless we average density matrices over time. The following theorem estab- 
lishes the uniqueness and convergence properties provided that the quantum walk 
satisfies a certain condition. Let us call a density matrix p strictly positive and write 
p > 0, if (x, px) = implies that x = 0. Next, let T be a linear operator acting 
on M.. We will call T strongly positive if for every density matrix p there exists an 
integer n > such that T n p > 0. 

(This definition is similar to a corresponding definition in the theory of Markov 
chains, in which it is shown that a stochastic matrix of a Markov chain is strongly 
positive if and only if the Markov chain is ergodic, that is, aperiodic and irre- 
ducible.) 

The multiplicity of an eigenvalue A is defined as dim ker (XI — T) . The rank of 
A is sup p>0 dim ker (XI — T) p . The eigenvalue is called simple if its rank equals 
1. 

Theorem 1.4. Let T be a strongly positive superoperator. Then, ( i) T has a simple 
eigenvalue 1. (ii) The corresponding eigenvector p s t is a strictly positive density 
matrix. (Hi) For every initial density matrix po, T n po —> p s t, as n — )■ 00. 

Proof is in Appendix. 

After the convergence to the stationary distribution is established, it is natural 
to ask for an estimate on the mixing time. First, let us define the mixing time for 
a non-unitary quantum walks. The definition is different from the definition for 
the unitary walks since no time-averaging is necessary. The measurement at time 
t finds the walking particle at the vertex v in state s with probability p t (x, p) = 
(x|T* (p) \x), where x denote the pair (v, s) and p is the initial density matrix. If T 
is strongly positive, then these probabilities converge to a limit p (x) = (x\p s t\x) , 
which does not depend on the initial density matrix. Hence, we can define the total 
variation distance as d (t, p) = J2 x eVxS \ { x \'^ t (p) \ x ) ~ P ( x )\- The correspond- 
ing mixing time can be defined as 

tmix (e) = sup inf [t : d (t', p) < e for all t' > i} . 

p 
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Unfortunately, while it is easy to see that the asymptotic behavior of T* is gov- 
erned by the spectral radius of T, it is difficult to estimate the mixing time because 
of the non-normality of operator T. The essential difficulty is that for such opera- 
tors it is hard to estimate the duration of the transient behavior. It is the same prob- 
lem that makes it difficult to estimate the mixing time for non-reversible Markov 
chains. 

(In one particular example of a non-unitary continuous-time walk on cycle this 
difficulty has been overcome and an estimate on the mixing time has been derived 



We consider several examples of unitary walks in this paper. The table sum- 
marizes results for unitary quantum walks on a complete graph, a cycle, and a 
hypercube. 



It appears from this table that the mixing time for quantum walks is of simi- 
lar order as that for the corresponding classical random walks. In particular, the 
unitary quantum walks do not allow a quadratic speedup over classical walks, in 
contrast to the results for the mixing time in JTJ. The reason for this difference is 
that the mixing time defined in (TJ restricts the initial distributions of the particle 
to the class of distributions concentrated on a particular vertex of a graph, while 
we allow for arbitrary initial distributions. Note that this result does not rule out 
that the quadratic speedup can be achieved by non-unitary quantum walks. Some 
evidence in favour of this conjecture can be found in [7 ] and |fT51 . 

The rest of the paper is organized as follows. In the next section, we apply 
bounds on mixing times to particular examples of quantum walks on the cycle, 
hypercube, and complete graph. The proofs of the theorems are relegated to Ap- 
pendix. 



Example 1. (Cycle) . 

Proposition 2.1. The mixing time for the Hadamard quantum walk on the n-cycle 
satisfies the following inequalities: 



in EH.) 



Complete graph with n vertices 
Cycle with n vertices 
Hypercube with 2 n vertices 




2. Examples 



—n 2 < t mix (e) < —n 2 log (n) , 



where c\ and C2 are positive constants. 
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Proof: The eigenvalues of the Hadamard walk on the cycle with n vertices were 
found in [1J. They are 

^-^ (-(**) ± V i+sin2 (£*))■ 

where k = 0, . . . , n— 1. In order to describe the eigenvectors, let Xk, < k < n—1, 
be functions in L 2 (Z n ) defined by the formula Xk = Sr=o ex P 

(2vri^) S r . Then 

all eigenvectors have the form v ® Xfe where i> is a 2- vector that depends on fc. 

Indeed, if S and S 1 * are the left and right shift operator on I? (Z n ), respectively, 
then we can write U as a 2-by-2 block matrix, with blocks U\\ and [/12 equal 
S/ y/2, and blocks U 2 i and J7 2 2 equal -S* /\/2 and S*/y/2, respectively. It follows 

that U (v ® Xfc) = ^4fcW ® Xfe, where 



Let 




Then, eigenvectors of ^ can be written as vf 1 ' = (c, 1) with c = 1— \/2il exp {—Imkjn) .The 
corresponding eigenvectors of C7 are ul^ (8) Xfc with eigenvalues a! = t[ for 
a = 1,2. 
Note that 

0(1,2) (.SM)\ 1 

p£ := arg I r^, I = ± arccos 

The smallest difference between occurs when A; = and 1 and it can be es- 
timated by c/n 2 for a suitable constant c. It follows that the relaxation time is 
t re i ~ cn 2 , and by Theorem ll. li the mixing time is 

tmix (e) < -^n 2 log {n) , 

with a certain constant C2 > 0. 

It is easy to estimate the overlap of eigenvectors that correspond to eigenvalues 
/3q and pV It is greater than 0.97 for all n. By Theorem [PI we have 

QED. 

Example 2. (Hypercube) 

The mixing time of the quantum walk on a hypercube was previously studied in 
lfT2l . and we use their setup in the definition of quantum walk. The quantum walk 



V2 u 
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on the hypercube is also analyzed in [5 ] with emphasis on the hitting time of the 
walk. 

Consider a hypercube graph (Z2)"' with 2™ vertices. We think about vertices as 
indexed by numbers from to 2™ — 1 in the binary representation with n digits. 
The edges of the graph are put between numbers that are different in one bit only. 
The set of states S consists of n elements. We consider the Grover quantum walk. 
That is, a particle at vertex v in state s goes to the vertex w which is different from 
vertex v only in the bit s. It remains in state s with amplitude 2/n — 1 and goes to 
state s' with amplitude 2/n. 

Proposition 2.2. The mixing time for the Grover quantum walk on the n-dimensional 
hypercube satisfies the following inequalities: 



Proof: The eigenvalues and eigenvectors of the Grover quantum walk on the n- 
dimensional hypercube were found by Moore and Russell in lfl2~1 . The eigenvalues 
are 



where k = 0, . . . , n. We describe eigenvectors below. For the convenience of the 
reader, we also give a short verification of the result . 

For each sequence t = (t±, £2, • • • , t n ) of and 1, define xt £ L 2 (Zg) Dv tne 

formula xt ■ ■ ■ ,x n ) = 2~ n l 2 (—\y^ tiXi . All eigenvectors of the matrix U 

(i) (i) 

have the form v t ® xt, where v t is an n- vector that depends on t, and j = 
l,...,n. 

Indeed, the unitary matrix U can be written as a n-by-n block matrix, in which 
the ij -th block is bSj if i 7^ j and aSj if i = j. Here a = 2/n — 1, b = 
2/n and : I? (Zg) — > L 2 (Z?j) is the shift operator which acts as follows: 
(Skf) (xi, • • • , x n ) = f (x\,... ,Xf. + 1, ... , x n ) , where addition is modulo 2. 
Note that S kX t = (-l)* fc Xt- 

A computations shows that U (v xt) = A. (v) ® xt, where A is an n-by-n 
matrix (which depends on t) with entries = [SijCt + (1 — dij) b] (— 1) i . In 
other form, A = D + bP, where D = -diag ((-l)' 1 , (-1)' 2 , . . . , (-1)*") , and 
P=|l,l,...,l)<(-l) tl ,(-l)* a ,...,(-l)*»|. 

It is easy to verify that the following vectors are eigenvectors of A. Let k be the 
number of non-zero entries in vector t. First, assume that n > k > 1 and define 
x = ±u/ and v r = x 1_ ' r . Then v = (v%, . . . , v n ) is an eigenvector of A with 



— n < t m ix (e) < — n\og (n) . 




2k 
n 



y/k(n — k) 



n 
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eigenvalue 

, 2k Jk(n - k) 
\f = 1 ± 2i-^ -. 

n n 

In addition, note that every non-zero vector v such that v r = if t r = and 
Ylr=i v r = 0' i s an eigenvector of A with eigenvalue 1. The set of such vectors 
form an eigenspace of dimension k — 1. Similarly, every non-zero v such that 
v r = if t T = 1 and Ylr=i *V = is an eigenvector of A with eigenvalue —1. The 
set of such vectors forms an eigenspace of dimension n — k — 1. 

For the case when t = 0, the vector (1, . . . , 1) is an eigenvector with eigenvalue 
1 and its orthogonal complement is eigenspace of —1. For t = (1, 1, 1) , the 
situation is reverse. 

By counting dimensions of eigenspaces, it is clear that these are all eigenvalues 
of matrix A. Since there are 2™ different choices of vector t, we also found all 
eigenvalues of matrix U. It follows that these eigenvalues are ±1, and for k = 

I, ...,n-l. 

From the formula for eigenvalues, the distance between distinct eigenvalues can 
be estimated from below as A > ^ . Hence, t re i < . 
By applying Theorem ll.il we find 

tmix {£) < — nlog(n) . 

For the lower bound, consider for simplicity the case of even n = 2m. (The case 
of odd n is similar.) Let v x ^ denote the value of function v G L 2 ((Z2)™ xZ„) 
on vertex x and state k, and consider the eigenvectors that correspond to eigenval- 
ues and X m+ i respectively. Then, it is easy to compute the overlap of these 
eigenvectors as y/l — (m + 1) / (2m 3 ) ~ 1 for large n. The distance between ar- 
guments of eigenvalues and A m+1 is approximately 2/n. Hence, by Theorem 

II. 2l we have the inequality 

QED. 

Example 3. (Complete graph) 

There are several ways to define a discrete-time walk on the complete graph 
with n vertices. We will consider the following variant. Let \S\ = n. Define the 
entries of the unitary matrix as follows. 

U (ws',vs) = 5 WS 8 vs i + (1 ~ bos') j ■ 

In words, let the particle start at vertex v in state s. Then at the next moment of time 
it will be at vertex w = s. The particle moves to state s' with amplitude (—2/n) 
for s' 7^ v. If s' = v, then the amplitude of the transition is (1 — 2/n) . 
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Proposition 2.3. The mixing time for the quantum walk on the complete graph 
satisfies inequalities: 



where c is a positive constant. 

Proof: We will show that the eigenvalues of U are 1, — 1, i and —i with multi- 
plicities n(n — l)/2, 1 + (n — l)(n — 2)/2, n — 1, and n — 1, respectively. Let 
X vs = ip ( v i s ) ■ Then, the action of U can be written as follows: 



where X T is the transposed matrix X, l n is the column n-by-1 vector that consists 
of all ones, and 1^ is the corresponding row vector. 

If X = X T and all columns of X sum to 0, then U(X) = X. This gives us an 
eigenspace of operator U with eigenvalue 1 and dimension n(n — l)/2. 

Similarly if X T = —X and all columns of X sum to 0, then U (X) = —X. This 
gives us an eigenspace of U with eigenvalue —1. In addition, U (I) = —I. Hence, 
the dimension of the eigenspace with eigenvalue — 1 is 1 + (n — l)(n — 2)/2. In 
order to find the eigenspaces with eigenvalues dbi, consider U 2 . It acts as follows: 



Let c\ , . . . , Cn and r\ , . . . , r n be arbitrary numbers satisfying the conditions J2 c i = 
rj = 0, and define Xij = r$ + Cj. Then U 2 (X) = —X. Hence, these matrices 
belong to the eigenspace of U 2 with eigenvalue —1. The dimension of this space is 
2n — 2. It follows that U has two eigenspaces of dimension n — 1 which correspond 
to eigenvalues i and —i, respectively. By counting dimensions we confirm that we 
have found all eigenvalues and eigenspaces of matrix U. 

The relaxation time is t re \ = l/y/2. By applying Theorems 11.21 and [TTT1 the 
mixing time satisfies inequalities 



for some positive c\ and C2 . 

Appendix A. Proofs of Theorems 

Proof of Theorem ll.lt Let V& be the eigenspace corresponding to eigenvalue 
Pk of operator U. Then, we can write 



_ ''mix — 






in 




fc=l 
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where ip k € V k , \\y>k\\ = 1, J2k \ c k\ 2 = 1- Then at time t, 

m 

ip(x,t) = ^2c k ip k (x) (A fc )*, 
k=\ 

and 

m 

k=l k+l 

Hence, 

j T-l m 

™ X] l^(»,*)| 2 <i* = ^\ck\ 2 Wk{x)\ 2 



T 

t=Q k=l 



+^J2c k ci^ k (x) <p t ( 



-. , x) (^-f - 1 



It follows that 

^ T-l m 

p(x) := lim -Vl^i.i)! 2 = V |c fc | 2 (x)| 2 , 

T->oo J z — ' L — ' 
t=l k=l 



and 



T 

xeG 



rp 

E— ( \- I \ ( X k X l) ~ 1 
C k Cl<p k (X) (fit (x) — p= — 

k+l XkX l ~ 1 



In order to bound this quantity, note that 

^2 \ C k~Wk {x)% (x)\ < Q(ipk,<fl) \Ck\ |Q| < |Cfc| |q| , 

xeG 



and, therefore, 



d(r)<-j; 



T ^ | e i(A-A) _ i| • 



Note that \ e i ^ k ~^ - l| > ^d((3 k ,(3i) where d(p k ,Pi) is the distance between 
/Sfc and /3; modulo 2-7T, that is, <i ((3 k ,/3i) := min s {\/3i — (3 k + 27rs|} . Let A denote 
mim^ ci {PkiPl) an d assume that < fi\ < fa < ■ ■ ■ < fi m < 2-7T. Then we can 
write d (Pk,Pl) > Ad (fe, /) , where d (k, I) = min {\k — l\ ,\k — I + m\ ,\k — I — m 
This inequality holds because the shortest arc of the circle between j3 k and fii con- 
tains d(k,l) non-overlapping intervals whose endpoints are fa and the length of 
each of these intervals is at least A. It follows that | eWk-Pi) - 1 1 > ^d (k, I) . 
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Next, 



< — y 



1 I I I 2 ■ I |2 

2 ( \ C k\ + |Q 



d(k,l) 



—y 



This sum can be estimated as follows: 

1 2 



l c fc| 



k^l 

We estimate 



d(M) 



< 



fc=i 



TA J) 



m / m - \ 

fc=l \Z=2 K ' '/ 



[m/2] 



Z=2 V ' ' fc=l 

and since ^Zfcli l c fc| 2 = 1> we nave 



< 2 log (2m) 



Hence, 



27T 

d(T) < ^ log (2m) 



W (e) < 2vrlog(2m) — . 



QED. 

Proof of Theorem ll.2t Let the initial function be <p = (ip + ip') / y/2 and let 

A := d (A, A') . Then we compute: 



d(T) = — 



1 - cos AT + 



sin AT sin A 



Q{W) 



1 — cos A 

Consider interval I r = A _1 7r(l/6 + 2r, 5/6 + 2r), where r is a non-negative 
integer. For every T G I r , sin (AT) > 1/2. Moreover, there is an integer T r E I r 
because \I r \ = A -1 (2/3) tt > 1. Note that the distance between T r and T r+ i is 
less than 3A _1 7r. 

For every integer T r we have 



d(T r ) > 



1 / sin A 
III I 1 - cos A 



> 



1 



4T r A 



(The second inequality holds because A < 2, and therefore sin A/ (1 — cos A) > 
A -1 .) 

It follows that if e > is sufficiently small then there exists an integer T r be- 
tween Q (ip, tp') /4e A - 3A -1 7r and Q (ip, if>') /4eA, such that d(T r ) > e. In par- 
ticular, if e < Q(^,VO/80,thenT r > QC0,VO /4eA— 3A -1 7r > Q(ip,tp') /8eA, 
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and we conclude 

t ■ (e) > » 

''mix \t) il g e /± 

QED. 

Proof of Theorem 11.41 The proof of (i) and (ii) is an application of results by 
Krein and Rutman from [9]. In this paper a cone K in a Banach space X is fixed 
and operator A G L (X) is called strongly positive if for every non-zero x € K, 
there is an integer n > 0, such that A n x is in the interior of K. Theorem 6.3 of 
this paper (on page 70 of the English translation) shows that if A is compact and 
strongly positive, then there exists one and only one eigenvector of A in the interior 
of K and the corresponding eigenvalue exceeds all others in absolute value. More- 
over, the proof of the theorem shows that this eigenvalue is simple. The claim of 
our theorem follows if we apply the Krein-Rutman theorem to the cone of positive- 
definite matrices. Indeed, by 11.31 there exists p st such that Tp s t = Pst- Since T 
is strongly positive, hence p s t is in the interior of K (i.e., strictly positive); by 
the Krein-Rutman theorem it is the only eigenvector in the interior of K, and its 
eigenvalue 1 is simple. 

For (iii), let Z be the space of Hermitian matrices with zero trace, Z = {p : tr (p) 
Then, TZ C Z and all the eigenvalues of T\ z are less than 1 in absolute value, 
because 1 is a simple eigenvalue and p st ^ Z. It follows that the spectral radius rj 
of T\z is smaller than 1. Hence, lim n __ +00 IKTIz)™!! 1 ^™ = V < 1, which implies 
that T n z — > for every z G Z. Since po — p st belongs to Z for every density 
matrix po, we conclude that T n po — > p s t for every pQ. QED. 
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